A condition for any realistic theory of quantum systems 
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In quantum physics, the density operator completely describes the state. Instead, in classical 
physics the mean value of every physical quantity is evaluated by means of a probability distribution. 
We study the possibility to describe pure quantum states and events with classical probability 
distributions and conditional probabilities and prove that the distributions can not be quadratic 
functions of the quantum state. Some examples are considered. Finally, we deal with the exponential 
complexity problem of quantum physics and introduce the concept of classical dimension for a 
quantum system. 



A peculiar aspect of the standard quantum formalism 
is interference, that is, every alternative is not associated 
with a probabiUty, but with a complex number. This 
characteristic is at the basis of some famous paradoxes 
as the Schrodinger cat [ij, and of the "sign problem" 
encountered in quantum Monte Carlo simulations 2i|, 
which gives an exponential growth of the numerical com- 
plexity. Since the birth of quantum mechanics, many 
physicists investigated the possibihty of its description 
in terms of classical probabilities. Indeed, a completely 
equivalent alternative to the Copenhagen interpretation 
is Bohni mechanics [3|, |j|, recently extended to quan- 
tum fields with a variable particle number Q. These 
approaches provide a realistic description of quantum 
systems. Bohm mechanics uses the same mathematical 
tools of the standard approach, as the wave-function, so 
it differs only in the interpretation and does not solve for 
example an important problem as the exponential com- 
plexity. The wave-function is not defined in the three- 
dimensional physical space, as a classical field, but in 
the representation space, thus the number of variables 
needed to define the classical state grows exponentially 
with the dimension of the physical system. Since quan- 
tum mechanics has a statistical interpretation, there is no 
a priori reason to exclude the possibility to describe the 
quantum states as statistical ensembles in a smaller space 
than the Hilbert space. In some cases, the Wiener func- 
tion provides such a dimensionality reduction [7| . It has 
some properties of a classical probability distribution in 
phase-space, but in general it can take negative values. 
When it is positive and only particular measurements 
are performed, a realistic statistical description in phase- 
space is possible, as in the case of continuous-variable 
teleportation experiments involving Gaussian states and 
quadrature measurements [iJ]. 

In order to circumvent the sign problem of the Wigner 
function, some alternative non-negative probability dis- 
tributions were introduced [8:, 9|, as the Husimi function. 
However, also in these cases, a realistic interpretation 
is not in general possible. Another example of phase- 
space distribution is the P-function, introduced by R. J. 
Glauber [l^. 

A realistic statistical description of quantum mechan- 
ics in a reduced phase-space is very important for its pos- 
sible consequences in Monte Carlo simulations of many- 



body systems. In the cases where it is possible to inter- 
pret the Wigner function as a probability distribution, as 
in the truncated Wigner approximation, the Monte Carlo 
approach allows to simulate efficiently the dynamics of 
atoms in degenerate bosonic gases [11[ and electrons in 
semiconductors [ij] • A general discussion of the method 
is also reported in Ref. flSj . 

In this Letter, we analyze the possibility to have a real- 
istic description for experiments involving general states 
and measurements and prove that a necessary condition 
has to be fulfilled. Using the non-conflicting hypothe- 
ses of the theorem, we finally introduce the concept of 
classical dimension of a quantum system. In this Letter 
we use the term "classical" as equivalent to "realistic", 
thus, a classical theory is a theory whose variables have 
definite values and for which it is possible to use the clas- 
sical rules of the probability theory. In this sense, Bohm 
mechanics is a classical theory and, consequently, every 
state is classical, i.e., has a realistic description. 

We associate each quantum state Jt/i) with a distribu- 
tion of probability p{X\'ip) in a suitable classical space, 
spanned by a set of variables X. This space can be more 
general than the phase-space and with higher dimension- 
ality. No a priori hypothesis on it is introduced. We want 
to give physical meaning to these variables, i.e., we as- 
sume that there is an underlying theory of the quantum 
system described by them. When the system is prepared 
to a state lip), X takes a value with a probability given 
by p. If a von Neumann measurement is performed, the 
state collapses to a state \(j>) with probability |(^|0)p. In 
terms of the realistic theory, we say that the system with 
coordinates X has a conditional probability P{(f>\X) to 
give the event </>. The probability of this event, given the 
state \tp), is obtained integrating P{(j)\X)p{X\ip) over X, 



M^)f^Tr[P^R^ 



dXPi^\X)p{X\^), (1) 



where P^ = \ip){ip\. The functions p{X\tp) and P{(I>\X) 
have to satisfy the conditions 



p{X\^) > 
dXp{X\ij) = 1, 



(2) 
(3) 



0<P{(f>\X) < 1. 



(4) 



for every [ip) and [0). 

A relation as Eq. ^ holds when we evaluate the prob- 
ability distribution of the position of a particle by means 
of the Wigner function. The probability to have a parti- 
cle in the spatial region ft is 



dqdpPin\q,p)W{q,p) 



where 



Pin\q,p) 



8{q-q)dq. 



(5) 



(6) 



Since < P{Vt\q^p) < 1, P it can be interpreted as a 
conditional probability. Thus, if the Wigner function is 
positive and only position measurements are involved, a 
realistic description in phase-space is possible |14| . 

Every probability distribution introduced so far is 
quadratic in the pure quantum state |'0), as the P, the 
Wigner and the Q functions (lO|. Indeed, the Wigner 
function was rigorously derived in Ref. [l5| using this 
property and other four assumptions. We will show that 
quadratic distributions can not be interpreted as proba- 
bility distributions of some realistic theory when general 
pure states and measurements are considered. More pre- 
cisely, by assuming that properties (HJH) are fulfilled for 
every ip and (16i | , we prove that the probability distribu- 
tions associated with pure quantum states are nonlinear 
functions of the density operator. This is the main result 
of our work. Note that Bohm mechanics is a perfectly 
consistent realistic hidden variable theory and the cor- 
responding probability distributions are not quadratic in 
the wave-function, as we will show. 

In general, a positive probability distribution linear in 
the density operator can be obtained by means of pos- 
itive operator- valued measurements (POVM) [l3|. We 
consider a pure state {ip) and define the associated prob- 
ability distribution with respect to a variable X as 



p{X\i;)=Tr[A{X)P^] 



(7) 



where P^ = \ip){'4'\ is a projector and A{X) is a generic 
Hermitian matrix which depends on X. X can be a vec- 
tor of continuous and/or discrete variables. p{X\'ip) must 
satisfy the properties (|2l3p . The first one is fulfilled if 
A{X) is positive definite, i.e. if its ei genv alues are non 
negative. The second one implies that [l3| 



/ 



dXA{X) = 1. 



(8) 



After these preliminary remarks, we can prove the theo- 
rem. Equations (fTI- l4|7|) are our starting hypotheses. We 
will show that they are conflicting. 
From Eqs. (|ll7p . we have 



Tr 



dXP{cb\X)A{X) ~ P^]P^ } = 0, 



for every \ip) and \<f>). Thus, 

p^= f dXP{(j)\x)A{x). 



(9) 



Since A{X) is positive definite and P{(j)\X) > 0, it is 
evident that 

P(0|X) ^0^ AiX) oc P^. (10) 

Thus, if P{(I)\X) ^ and P((/)'|X) ^ 0, then = <j>' . We 
can define the following one- valued function 

X-X^ X{X) such that P{x{X)\X) ^ 0. (11) 

The function x(^) spans the entire Hilbert space, but it 
is not necessarily invertible. It is possible to introduce 
an auxiliary set Y{X) of variables in order to have the 
bijective mapping X ^^ (x, i^). Condition (|10p implies 
that 



P(0|x,n = oif0^x 



(12) 



A{x,Y) = a{x,Y)P^, (13) 

where a(x, Y) > 0. Equations (|9ll2ll3p give 



VxdYP{cP\x,Y)a{(b,Y) = l, 



(14) 



for every (p. Since P{(p\x,Y) is different from zero in a 
set with zero measure [Eq. (IT^ ] and it is smaller or equal 
to 1, J dYa{(p,Y) must be infinite for every </>, but this 
is impossible, because of the normalization condition ([5]). 
D 

Since the properties ^i^ are the minimal necessary 
conditions for a realistic theory, we have to discard the 
linear assumption, i.e., Eq. ([7]). This is sufficient to re- 
move the contradiction, as shown below, when we will 
explicitly write the probability density and conditional 
probability for the Bohm theory. 

Now, we consider some examples to illustrate our re- 
sult. The P, Wigner and Q functions can be defined 
for systems described by boson creation and annihila- 
tion operators, a^ and a, respectively [lO[. It well-known 
that the P and Wigner function cannot be probability 
distributions. The first one is highly singular for par- 
ticular states, as squeezed states or superposition of co- 
herent states. The second one is well-defined for ev- 
ery state, but can assume negative values, i.e. equa- 
tion ^ is not satisfied. The Q-function requires a more 
detailed discussion. It is smooth, positive and normal- 
ized, thus our theorem says that equation ^ cannot 
be satisfied. Consider a one-mode system and denote 
the Q-function corresponding to a state \ip) by Q{a\ip), 
where a is a complex number. The expectation value 
of an observable M{d\d) in anti-normal form is given 
by the average of the function M{a* , a) with respect to 
Q{a\'4>), The expectation value of the number of particles 



is (a'^a) = (aa'l') - 1 = J cP a{\a\'^ - l)Q{a\'ip), thus we 
have that "^"^^q nP {n\a) — jap — 1, where P{n\a) is the 
conditional probability of measuring n particles for the 
phase space state a. This equation implies that the con- 
ditional probability has to be negative for some values of 
n and |ap < 1, in agreement with our theorem. 

Now, we consider the case of a two dimensional Hilbert 
space. It is possible to construct a probability distribu- 
tion pinltp) with the following 3 matrices [l8|: Ai = 



^(Ti and A3 



1 



1^3 



-^(Ti, where a^ are the Pauli matrices. The distribu- 
tion has three values, is positive and normalized. Also 
in this case, the conditional probability P{(j)\n) cannot 
satisfy Eq. ((4]). For example, when cj) = (1,0), Equa- 
tion HI) is fulfilled if and only if P((?!)|l) = 3/2 and 
P(0|2) = P((/.|3) = 0. 

We can define another probability distribution 
as p{e,^\^) = ^Tr[\e,^}{e,(l)\p], where |0,0) = 
cos(6l/2)|0)-|-e''^sin(6l/2)|l) and |0}, |1) are two orthonor- 
mal vectors. 6 and are the polar coordinate of a point 
of the Bloch sphere. In the integrals with respect to these 
coordinates, we use the measure svaddddcj). The distribu- 
tion is positive and normalized. For jV'} = |0), we have 
p(6»,0|O) = [cos(6'/2)]2/27r. The probability of obtaining 
the state |0) is equal to 1, thus 



27r 



d(j) I sin 6'P(0|6l, 0)p(6l, 0|O) = 1 (15) 



We assume ab absurdo that P(O|0, 0) is positive, then 
the equality is satisfied only if P(0|6', 0) = 1, for 9 ^ -k. 
But this implies that the probability to obtain |0) when 
the system is in the state |1) is different from zero, since 
the corresponding probability distribution p{0^(j)\\) — 
[sin(0/2)]^/27r is zero only for 9 = Q. This is absurd 
because |0) and |1) are orthogonal. All the previous ex- 
amples show that at least one property required by the 
theorem is not satisfied. 

Discarding the linear hypothesis, we investigate how to 
construct a probability distribution associated with ev- 
ery density operator which satisfies the conditions ([TJIl]). 
The Bohm mechanics provides a simple example of true 
probability. For sake of simplicity, here we consider only 
the case of a fixed number of particles [sl, 14] . Recently, 
the theory has been extended to quantum fields and ac- 
counts for creation and annihilation of particle [5|. The 
dynamical variables are a multi-particle wave-function x 
and the coordinates x in the configuration space. If the 
quantum system is in a pure state ip, the variable x is 
merely equal to ip and the distribution of the coordinates 
X is \'4){x)\'^ , thus 



p(^,xl'0) = IV'(^)l^'5(x-V'), 



(16) 



where V' is the state in the Hilbert space of the quan- 
tum system. This function is obviously non-quadratic 
in ip^ because of the Dirac delta. In Bohm mechanics, 
the position measurement gives the variables x as result. 



thus, the conditional probability of finding the system in 
a volume f7 in the representation space is 



P{^\x,x) 



dxQd{xQ — x). 



(17) 



Other measurements, for example of momentum, can be 
performed by a suitable unitary evolution and a subse- 
quent measurement of position. During the evolution, 
the coordinates a' go to a new value which is function of 
X and ij}. This enables to obtain a positive conditional 
probability for every measurement from Eq. p7p . 

Bohm mechanics is a free-dispersion theory, i.e., the 
hidden variables fix exactly the results of measurements. 
It corresponds to have conditional probabilities equal to 
or 1. This characteristic is not necessarily required by 
a realistic theory. Indeed, the minimal requirements are 
Eqs. imi). Discarding the free-dispersion hypothesis, the 
simplest probability distribution for a pure state ij) is 



p(xlV') =='5(x-V'), 



(18) 



where x is the variable of the classical system. The 
corresponding conditional probability for the event </> is 

p{<p\x) = mx)?- 

Although these two examples sound trivial, they show 
that properties ([T][4|) are not conflicting. At this point, 
we raise a question regarding the exponential complex- 
ity of quantum mechanics. In the case of Eq. (US]), the 
dimension of the phase space spanned by the variables 
X is obviously the dimension of the Hilbert space of '0- 
It is well-known that the dimension of the Hilbert space 
grows exponentially with the physical dimension of the 
system. For example, it is 2^ for N spins 1/2. Thus, also 
the number of variables x of the classical theory has an 
exponential growth. Every known hidden variable the- 
ory has this feature. However, if we discard some re- 
quired property, as positivity, the dimension of the X 
space can be considerably reduced. For example, the 
Hilbert space dimension of one boson mode is infinite, 
since the number of particles goes from zero to infinity, 
but the corresponding Wigner functions have only two 
variables, although we cannot regard these variables as 
describing a classical system, because of the negativity of 
the Wigner function. In general, for a Hilbert space with 
dimension M, it is possible to find quasi-probability dis- 
tributions where X can assume only M x M values, that 
is, if the Hilbert space dimension is finite, we can have 
a space of X with dimension equal to zero. So, on one 
hand we have an example of true probability distribu- 
tion on a space with the same dimension of the Hilbert 
space, on the other one, we can have quasi-probability 
distributions on a space with a considerably lower di- 
mension with respect to the Hilbert space. The question 
we raise is the following: if the dimension of the Hilbert 
space is D quant , which is the lowest dimension Ddass for 
X in order to fulfiU the conditions dTJH])? Equation p^ 
says that Ddass ^ Dq^ant- The introduction of a prob- 
ability distribution p{X\'ip) and a conditional probabil- 
ity P{(f)\X) for each event allows to place the question of 



the exponential complexity nature of quantum mechanics 
onto a clear and well-defined ground and the proof that 

Dclass = Dquant Or that Ddass < Dguant WOuld Cast UCW 

light upon this question. The evaluation of Ddass is not 
a trivial problem and there is no evident reason for as- 
suming it equal to Dquant- We have demonstrated that 
the solution of the problem requires to discard probabil- 
ity distributions quadratic in the wave- function. Note 
that we have not dealt with dynamical considerations, 
but only with the problem of a classical representation of 
quantum states. 

In conclusion, we have studied the possibility of a clas- 
sical description of quantum states and events by means 
of probability distributions and conditional probabilities, 
that must satisfy the properties ([TJIl]). We have demon- 
strated that the probability distribution for a pure state 
has to be a nonlinear function of the density operator. 
We have illustrated the proof with some examples, as the 
P, Wigner and Q functions, and two cases of POVM for 
a two-state system. We have shown that for Bohm me- 



chanics the distributions of probability are not quadratic 
functions of the pure quantum state, as required by our 
theorem to any realistic theory. Finally, the concept of 
classical dimension Ddass of a quantum system has been 
introduced. Every known hidden variable theory has not 
a phase-space dimension smaller than the Hilbert space 
dimension Dquant- We conclude with a non-trivial ques- 
tion. Does a classical theory exist whose phase space 
dimension is much smaller than D, 



Dr 



<^D, 



quant 



quant'^ that is, is 
? Bohm mechanics could be considered 



as a pure interpretation curiosity, since it uses essentially 
the same tools of quantum mechanics, whereas a theory 
with much lower dimensionality would have important 
implications. We have put this problem onto a well- 
defined ground and its solution could clarify the nature 
of the exponential complexity of quantum mechanics. 
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